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ABSTRACT 

The  authors  conduct  an  analysis  of  the  dynamics  of  secondary  eyewall  formation  in  two  modeling  frame¬ 
works  to  obtain  a  more  complete  understanding  of  the  phenomenon.  The  first  is  a  full-physics,  three- 
dimensional  mesoscale  model  in  which  the  authors  examine  an  idealized  hurricane  simulation  that 
undergoes  a  canonical  eyewall  replacement  cycle.  Analysis  of  the  mesoscale  simulation  shows  that  secondary 
eyewall  formation  occurs  in  a  conditionally  unstable  environment,  questioning  the  applicability  of  moist- 
neutral  viewpoints  and  related  mathematical  formulations  thereto  for  studying  this  process  of  tropical  cyclone 
intensity  change.  The  analysis  offers  also  new  evidence  in  support  of  a  recent  hypothesis  that  secondary 
eyewalls  form  via  a  progressive  boundary  layer  control  of  the  vortex  dynamics  in  response  to  a  radial 
broadening  of  the  tangential  wind  field. 

The  second  analysis  framework  is  an  axisymmetric,  nonlinear,  time-dependent,  slab  boundary  layer  model 
with  radial  diffusion.  When  this  boundary  layer  model  is  forced  with  the  aforementioned  mesoscale  model’s 
radial  profile  of  pressure  at  the  top  of  the  boundary  layer,  it  generates  a  secondary  tangential  wind  maximum 
consistent  with  that  from  the  full-physics,  mesoscale  simulation.  These  findings  demonstrate  that  the 
boundary  layer  dynamics  alone  are  capable  of  developing  secondary  wind  maxima  without  prescribed  sec¬ 
ondary  heat  sources  and/or  invocation  of  special  inertial  stability  properties  of  the  swirling  flow  either  within 
or  above  the  boundary  layer.  Finally,  the  time-dependent  slab  model  reveals  that  the  simulated  secondary 
wind  maximum  contracts  inward,  as  secondary  eyewalls  do  in  mesoscale  models  and  in  nature,  pointing  to 
a  hitherto  unrecognized  role  of  unbalanced  dynamics  in  the  eyewall  replacement  cycle. 


1.  Introduction 

Secondary  eyewall  formation  (SEF)  is  widely  recog¬ 
nized  as  an  important  research  problem  in  the  dynamics 
of  mature  tropical  cyclones,  but  as  of  yet  there  is  not 
a  consensus  on  the  phenomenon’s  fundamental  physics. 
The  conceptual  and  empirical  linkage  between  secondary 
eyewalls  (SEs)  to  hurricane  intensity  change  (Willoughby 
et  al.  1982,  1984;  Kuo  et  al.  2008)  and  storm  growth 
(Maclay  et  al.  2008)  has  fostered  renewed  interest  in 
SEF  (Houze  et  al.  2007),  with  substantial  efforts  cur¬ 
rently  under  way  in  developing  statistical  forecast¬ 
ing  tools  for  the  operational  community  (Kossin  and 
Sitkowski  2009).  However,  as  of  today,  such  forecasting 
instruments  tend  to  rely  more  on  empirical  relations 
than  in  the  understanding  of  the  physical  processes  that 
lead  to  SEF. 
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Previous  studies  of  the  physics  of  SEF  have  sum¬ 
moned  a  variety  of  physical  processes.  These  include  the 
interaction  of  the  tropical  cyclone  core  with  nearby 
weaker  vorticity  anomalies  (Kuo  et  al.  2004,  2008);  a 
finite-amplitude  wind-induced  surface  heat  exchange 
(WISHE;  Emanuel  1994,  1995,  1997)  instability  trig¬ 
gered  by  external  forcing,  such  as  with  the  interaction  of 
a  tropical  cyclone  with  a  baroclinic  eddy  in  the  upper 
troposphere  (e.g.,  Nong  and  Emanuel  2003);  the  hy¬ 
pothesized  impact  of  vortex  Rossby  waves  (VRWs) 
on  the  tropical  cyclone  inner  core  (Montgomery  and 
Kallenbach  1997;  Qiu  et  al.  2010;  Martinez  et  al.  2010; 
Menelaou  et  al.  2012;  Abarca  and  Corbosiero  2011), 
or  anisotropic  upscale  energy  cascades  (Terwey  and 
Montgomery  2008,  hereafter  TM08).  It  has  been  pro¬ 
posed  recently  that  SEF  results  from  high  generation  and 
accumulation  of  potential  vorticity  and  that  VRWs  are 
not  a  contributing  factor  for  SEF  (Judt  and  Chen  2010). 
Another  proposal  is  that  the  sustained  latent  heating 
outside  the  primary  eyewall  can  lead  to  SEF  (Rozoff  et  al. 
2012)  in  association  with  a  radial  expansion  of  the  kinetic 
energy  and  an  increase  of  the  axisymmetric  efficiency 
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there  as  defined  by  Schubert  and  Hack  (1982)  and  Hack 
and  Schubert  (1986). 

Only  recently  has  a  clear  dynamical  link  been  made 
between  the  overarching  mechanisms  of  tropical  cyclone 
intensification  and  the  physics  of  SEF  (Huang  et  al.  2012). 
This  realization  followed  soon  after  the  development  of 
the  paradigm-shifting  study  of  Smith  et  al.  (2009),  in 
which  it  was  shown  that,  from  the  perspective  of  the  az- 
imuthally  averaged  system-scale  flow,  the  spinup  of  the 
swirling  circulation  occurs  via  two  mechanisms. 

The  first  spinup  mechanism  occurs  above  the  fric¬ 
tional  boundary  layer  by  the  conventional  picture  of 
inflowing  rings  of  fluid  that  (approximately)  materially 
conserve  their  azimuthally  averaged  absolute  angular 
momentum  M  =  rv  +  (1/2) /r2,  where  r  is  the  radius  from 
the  circulation  center,  v  is  the  azimuthally  averaged 
tangential  velocity,  and /is  the  Coriolis  parameter.  The 
inflowing  air  occurs  primarily  below  the  midtroposphere 
and  results  from  the  radial  gradient  of  buoyancy  gen¬ 
eration  rate  associated  with  the  aggregate  positive  latent 
heating  in  the  deep  convective  clouds  that  span  the 
troposphere.  The  latent  heating  rate  of  these  clouds  is 
maximized  near  the  midtroposphere  (~7km).  Since  the 
frictional  torque  on  fluid  rings  above  the  boundary  layer 
is  negligible,  and  assuming  pressure  torques  from  hori¬ 
zontal  eddy  motions  are  negligible,  M  will  be  materially 
conserved  there.  Then,  the  foregoing  relationship  be¬ 
tween  M  and  v  implies  that  as  deep  convection  draws  in  rings 
of  fluid,  the  tangential  velocity  will  increase.  This  mechanism 
has  been  articulated  previously  by  many  authors  (e.g., 
Ooyama  1969;  Carrier  1971a;  Willoughby  1979;  Shapiro 
and  Willoughby  1982;  Willoughby  1995).  Using  a  non¬ 
hydrostatic,  convection-representing  model,  recent  work 
has  affirmed  that  this  mechanism  may  be  interpreted  in 
terms  of  axisymmetric  balance  dynamics  and  has  used  this 
mechanism  to  explain  why  the  vortex  expands  in  size  (Smith 
et  al.  2009;  Bui  et  al.  2009;  Fudeyasu  and  Wang  2011). 

The  second  spinup  mechanism  occurs  within  the 
frictional  boundary  layer,1  when  the  depletion  of  M  by 
nonconservative  forces  is  surpassed  by  the  inward  radial 
advection  of  M  associated  with  the  strong  frictional  in¬ 
flow.  Although  M  is  not  materially  conserved  in  the 
boundary  layer,  positive  tangential  wind  tendency  (resulting 


1In  this  paper  we  follow  Smith  at  al.  (2009)  and  Smith  and 
Montgomery  (2010)  and  adopt  a  dynamical  definition  of  the 
boundary  layer,  using  the  term  “boundary  layer”  to  describe  the 
shallow  layer  of  relatively  strong  inflow  near  the  sea  surface  and 
that  arises  largely  because  of  the  frictional  disruption  of  gradient 
wind  balance  near  the  surface.  It  will  be  shown  below  for  the  nu¬ 
merical  simulation  analyzed  in  detail  here  that  the  boundary  layer 
so  defined  is  about  1  km  deep  in  the  primary  eyewall  region  of  the 
hurricane  and  somewhat  deeper  (—1.5  km)  in  the  region  of  SEF. 


in  large  tangential  wind  velocity)  can  be  achieved  there  if  the 
radial  inflow  is  sufficiently  large  to  bring  the  air  parcels 
to  small  radii  with  minimal  loss  of  M.  This  spinup  mecha¬ 
nism,  while  coupled  to  the  interior  flow  via  the  radial 
pressure  gradient  at  the  top  of  the  boundary  layer,  is  tied 
fundamentally  to  the  dynamics  of  the  boundary  layer, 
where  the  flow  is  not  in  gradient  wind  balance  over  a 
substantial  radial  span  (Smith  et  al.  2009;  Smith  and 
Montgomery  2010).  Whereas  this  second  mechanism  cannot 
be  captured  by  axisymmetric  balance  dynamics,  numerical 
model  studies  and  recent  observations  show  it  to  be  the  es¬ 
sential  pathway  for  the  generation  of  the  maximum  tan¬ 
gential  winds  on  the  cyclone  scale  (Zhang  et  al.  2001;  Smith 
et  al.  2009;  Smith  and  Thomsen  2010;  Sanger  et  al.  2013).2 

Motivated  by  these  developments,  using  the  model- 
observation-consistent  dataset  of  Wu  et  al.  (2012),  Huang 
et  al.  (2012)  proposed  that  the  boundary  layer  force 
imbalance  in  the  radial  direction  and  its  radial  structure 
play  in  an  important  role  in  initiating  and  sustaining 
convective  activity  radially  outward  of  the  primary 
eyewall.  They  were  the  first  to  conceptualize  and  demon¬ 
strate  SEF  as  the  culmination  of  a  sequence  of  structural 
changes  in  a  mature  hurricane.  According  to  the  new 
model,  the  sequence  begins  with  a  broadening  of  the  tan¬ 
gential  wind  field  (via  the  first  spinup  mechanism)  followed 
by  an  increase  of  boundary  layer  inflow  outside  of  the 
primary  eyewall  and  the  development  of  supergradient 
winds  in  this  region  near  the  top  of  the  boundary  layer.  The 
latter  is  associated  with  horizontal  convergence  and  the 
emergence  of  the  deep  convection  that  form  the  roots  of 
the  SE.  They  suggested  also  that  simple  diagnostics  in¬ 
volving  the  agradient  force  could  prove  helpful  in  opera¬ 
tional  practice  to  help  diagnose  SEF. 

In  the  context  of  the  new  SEF  model,  we  believe  it  is 
useful  to  further  analyze  SEF  in  the  idealized  full-physics, 
three-dimensional  simulation  generated  by  TM08.  Al¬ 
though  the  TM08  simulation  was  used  originally  to  ex¬ 
plore  the  hypothesis  of  anisotropic  inverse  cascades 
mentioned  previously,  this  same  simulation  will  be  used 
here  to  test  the  new  model  advanced  by  Huang  et  al. 
(2012).  Since  this  simulation  has  been  further  demon¬ 
strated  to  furnish  a  canonical  SE  (Terwey  et  al.  2013), 
the  simulation  is  believed  to  be  an  ideal  laboratory 
for  such  a  test.  Because  boundary  layer  dynamics  are 
hypothesized  to  be  a  central  ingredient  in  both  the  es¬ 
tablishment  of  convective  activity  and  the  emergence  of 
a  secondary  wind  maximum,  the  TM08  simulation  will 
be  used  in  conjunction  with  a  simple  slab  boundary  layer 


2  A  review  of  the  most  prominent  tropical  cyclone  intensification 
paradigms  over  the  past  five  decades  is  presented  by  Montgomery  and 
Smith  (2013). 
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model  to  isolate  the  pure  boundary  layer  response  to  an 
imposed  wind  field  at  the  top  of  the  boundary  layer. 

An  outline  of  the  remaining  parts  of  the  paper  is  as 
follows:  Section  2  describes  the  two  numerical  models 
used  in  this  study  and  summarizes  the  analysis  meth¬ 
odology  employed.  Section  3  reviews  the  evidence  of  the 
canonical  SE  in  the  TM08  simulation.  Section  4  presents 
kinematic  characteristics  of  the  SEF  event.  Section  5 
presents  a  dynamical  analysis  of  the  SEF  event.  Section 

6  uses  the  tangential  wind  at  the  top  of  the  simulated 
boundary  to  compare  the  predictions  between  the  slab 
boundary  layer  model  and  the  TM08  simulation.  Section 

7  presents  our  conclusions. 


2.  Models  and  methodology 

a.  The  mesoscale  model 

This  paper  analyzes  the  SEF  simulation  initially 
studied  by  TM08  and  recently  revisited  by  Terwey  et  al. 
(2013)  that  used  the  Regional  Atmospheric  Modeling 
System  (RAMS;  Pielke  et  al.  1992;  Cotton  et  al.  2003). 
Briefly,  RAMS  is  a  three-dimensional,  nonhydrostatic 
numerical  modeling  system  with  two-way  interactive 
nesting  capabilities  (Clark  and  Farley  1984).  The  model 
simulation  uses  a  surface  flux  parameterization  for  mo¬ 
mentum  and  heat  based  on  the  Louis  (1979)  scheme;  the 
radiation  scheme,  which  includes  cloud  microphysical 
interactions,  introduced  by  Harrington  (1997);  the  subgrid- 
scale  turbulence  scheme,  based  on  Smagorinsky  (1963) 
and  that  incorporates  the  modifications  from  Lilly  (1962) 
and  Hill  (1974);  and  the  microphysical  scheme,  with  seven 
species  of  hydrometeors  by  Walko  et  al.  (1995). 

The  RAMS  simulation  is  an  idealized  integration  on 
an /plane  (centered  at  15°N),  with  constant  sea  surface 
temperature  of  28°C.  The  simulation  is  configured  with 
30  vertical  levels  with  vertical  grid  spacing  varying 
gradually  from  300  m  near  the  surface  to  1800  m  near  the 
top  of  the  domain  (at  about  26-km  altitude).  The  hori¬ 
zontal  grid  structure  is  nested  on  three  domains,  with 
horizontal  grid  spacings  of  24,  6,  and  2  km  and  168, 170, 
and  251  grid  points,  respectively.  The  respective  domain 
sizes  are  4032,  1020,  and  502  km.  Each  inner  grid  is  lo¬ 
cated  at  the  center  of  the  parent  grid.  The  outermost  grid 
is  cyclic  in  the  horizontal  direction  (to  help  constrain 
mass  conservation  in  the  domain).  Given  the  size  of  the 
largest  domain  and  the  relatively  small  size  of  the  sim¬ 
ulated  vortex,  it  is  thought  that  there  is  minimal  in¬ 
teraction  of  the  storm  with  itself  (TM08). 

The  simulation  is  initialized  with  a  quiescent  envi¬ 
ronment  in  which  the  vertical  thermodynamic  structure  is 
given  by  Jordan’s  (1958)  summertime  Caribbean  sound¬ 
ing.  The  initial  vortex  employed  is  a  weak  (10ms_1) 


mesoscale  cyclonic  vortex  in  gradient  and  hydrostatic 
balance,  located  at  the  center  of  the  three  domains.  The 
initial  vortex  reaches  its  maximum  tangential  wind  ve¬ 
locity  at  4-km  height  and  has  a  radius  of  maximum  winds 
(RMW)  of  75  km.  The  initial  conditions  include  also 
a  positive  moisture  anomaly  with  a  mono  tonic  increase  in 
water  vapor  mixing  ratio  by  1.3  g  kg-1  from  the  envi¬ 
ronment  to  the  center  of  the  initial  vortex.  Further  details 
of  the  RAMS  simulation  and  the  experimental  setup  can 
be  found  in  TM08  and  in  Montgomery  et  al.  (2006). 

Throughout  this  study  we  present  azimuthal  averages 
of  the  simulated  variables.  The  averages  are  obtained 
after  interpolating  the  RAMS  data  from  its  original 
Arakawa  C  grid  in  Cartesian  coordinates  into  a  cylindri¬ 
cal  grid,  whose  center  is  determined  as  the  centroid  of  the 
potential  vorticity  field  in  the  lowest  7.3  km  of  the  do¬ 
main.  All  of  the  results  presented  for  specific  times  cor¬ 
respond  to  2-h  averages,  centered  at  the  noted  time.  For 
the  present  analysis,  temporal  averages  are  preferred 
over  instantaneous  fields  because  the  time-averaging 
procedure  reduces  the  effects  of  transient  flow  features, 
such  as  inertia-gravity  waves  and  individual  convective 
cells  that  are  common  features  in  the  analyzed  dataset. 


b.  The  time- dependent  boundary  layer  model 

To  address  the  question  of  whether  boundary  layers 
can  generate  secondary  wind  maxima  in  response  to  an 
imposed  hurricane  wind  field  above  without  secondary 
wind  maxima,  a  slab  boundary  layer  model,  based  on 
that  described  in  Shapiro  (1983),  is  integrated  in  radius 
and  time.  The  model  solves  the  following  cylindrical 
radial  and  tangential  momentum  equations: 


d(u) 

~df 


d(u)  ^  (v)2 


CT 


~(u)^y  +  ^ -  +  f(v }  -  -f  (u)((uY  +  (vY) 


2x1/2 


(i) 


(2) 


where  u  and  v  are  the  radial  and  tangential  wind  ve¬ 
locities,  z  is  the  vertical  coordinate,  t  is  time,  p  is  the  air 
density,  p  is  pressure,  CD  is  the  drag  coefficient,  h  is 
the  prescribed  constant  boundary  layer  height,  V2  is  the 
horizontal  Laplacian  operator  in  cylindrical  coordinates, 
K  is  the  horizontal  eddy  diffusivity,  angle  brackets  denote 
azimuthal  average,  and  the  rest  of  the  symbols  have  their 
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standard  meanings.  In  Eqs.  (1)  and  (2),  the  terms  on  the 
left-hand  side  are  the  local  time  rate  of  change  of  the 
radial  and  tangential  velocity  components,  and  those  on 
the  right-hand  side  are,  from  left  to  right,  respectively, 
radial  advection,  centrifugal,  Coriolis,  frictional,  and  ra¬ 
dial  diffusion  forces.  The  radial  momentum  equation  [Eq. 
(1)]  contains  also  the  radial  pressure  gradient  force, 
which  is  the  driving  force  for  the  boundary  layer  flow. 
For  simplicity,  the  slab  model  neglects  vertical  advection 
(Shapiro  1983). 

As  is  standard  in  boundary  layer  calculations  (e.g., 
Smith  1968;  Carrier  1971b),  the  tangential  wind  at 
boundary  layer  top  is  used  to  force  the  slab  boundary 
layer  model  through  the  radial  pressure  gradient  term  in 
Eq.  (1).  The  tangential  wind  at  boundary  layer  top  is  assumed 
to  be  in  gradient  wind  balance  with  the  pressure  field. 

To  complete  the  specification  of  the  slab  model  a  suit¬ 
able  formulation  of  the  surface  drag  coefficient  is  required. 
After  the  observational  study  of  Bell  et  al.  (2012a),  based 
on  the  Coupled  Boundary-Layer  Air-Sea  Transfer  ex¬ 
periment  (CBLAST;  Black  et  al.  2007;  Drennan  et  al. 
2007;  French  et  al.  2007;  Zhang  et  al.  2008),  we  compute 
the  surface  drag  coefficient  following  Thomsen  et  al. 
(2013)  as  CD  =  0.7  X  10^3  +  1.4  X  10“3  [1  -  exp(- 
0.055 |u|)].3  For  the  integrations  presented  here,  a  con¬ 
stant  radial  eddy  diffusivity  of  K  =  1  X  104m2s_1  is 
used.  This  relatively  large  value  is  chosen  as  a  crude,  but 
plausible,  parameterization  of  the  shear  instability  and 
mixing  processes  at  the  inner  edge  of  the  eyewall.  This 
diffusivity  serves  also  to  regularize  the  numerical 
boundary  layer  solution  that,  in  the  absence  of  radial 
momentum  diffusion,  tends  to  produce  discontinuous 
“shocks”  in  both  the  radial  and  tangential  velocity  fields 
near  and  just  outside  the  RMW  (Shapiro  1983). 

The  magnitude  of  the  chosen  radial  eddy  diffusivity  is 
not  unlike  the  values  estimated  from  in  situ  observations 
by  Zhang  and  Montgomery  (2012),  who  reported  values 
as  high  as  K  =  3.43  X  103m2s_1  (their  Table  1)  in 
Hurricane  David  (1979).  The  magnitude  of  the  radial 
eddy  diffusivity  is  comparable  also  to  that  used  by 
Rotunno  and  Emanuel  (1987,  their  section  3a)  in  their 
axisymmetric  modeling  study  (attaining  values  as  large 
as  1  X  104m2s_1).  Qualitatively  similar  slab  boundary 
results  to  that  described  below  are  obtained  with  de¬ 
creased  eddy  diffusivities  as  low  as  6  X  103m2s_1. 

The  model  is  integrated  in  time  and  radius,  employing 
a  constant  Coriolis  parameter  evaluated  at  15°N  to 


3  Sensitivity  analysis  using  the  Deacon’s  formula  (Roll  1965)  and 

a  constant  CD  =  2.0  X  10-3  demonstrate  robustness  of  the  results, 
with  strong  qualitative  and  quantitative  agreement  among  a  cor¬ 
responding  suite  of  integrations  (not  shown). 


match  the  value  used  in  the  RAMS  simulation.  A 
constant  boundary  layer  depth  of  2  km  is  used  to  be 
broadly  consistent  with  the  boundary  layer  depth  in  the 
region  of  SEF  diagnosed  from  the  RAMS  simulation 
(see  upcoming  Figs.  2a-c  and  discussion).  With  the 
given  choices  of  CD  and  h,  the  vertically  averaged  slab 
boundary  layer  drag  force  magnitude  in  the  tangential 
direction  (between  75  and  150  km  of  radius)  is  2.41  X 
10-3ms-2.  In  comparison,  the  vertically  averaged 
tangential  drag  force  in  the  RAMS  simulation  (calcu¬ 
lated  as  the  residual  of  the  average  forces  in  the  lowest 
1138  m,  between  75-  and  150-km  radius,  at  hour  23 — 
during  SEF,  as  shown  later)  has  a  magnitude  of  1.38  X 
10-3ms-2,  with  local  values  reaching  magnitudes  of 
5.00  X  10  3  m  s  2.  In  other  words,  the  magnitude  of  the 
tangential  frictional  force,  as  well  as  the  rest  of  the 
parameter  choices  in  the  slab  boundary  layer  model, 
falls  within  the  range  of  values  obtained  in  the  RAMS 
simulation. 

3.  Secondary  eyewall  in  the  RAMS  simulation 

To  set  the  scene  for  the  upcoming  diagnostic  analyses, 
we  review  first  the  main  evidence  of  SEF  in  the  TM08 
simulation.  Figure  2  of  TM08  presents  a  time  series  of 
minimum  pressure,  RMW  (at  the  lowest  model  level  of 
147  m),  and  maximum  azimuthal-mean  tangential  winds 
(at  1  km)  for  the  220-h  numerical  simulation.  The  sim¬ 
ulated  storm  evolves  from  a  weak  mesoscale  convective 
vortex  to  a  mature  hurricane  that  subsequently  un¬ 
dergoes  a  canonical  eyewall  replacement  cycle.  In  the 
figure,  SEF  is  characterized  by  a  sudden  expansion  of  the 
RMW  (from  radius  35  to  80  km  right  before  hour  180), 
which  is  associated  with  a  weakening  of  the  tangential 
wind  in  the  primary  eyewall  and  an  intensification  of  the 
tangential  wind  in  the  SEF  radial  region.  Following 
TM08,  we  denote  hour  156  as  hour  0  to  more  easily  de¬ 
scribe  the  kinematics  and  dynamics  of  SEF. 

Figure  1  shows  radius-time  diagrams  of  the  azimuthal- 
mean  tangential  and  vertical  velocity  at  a  height  of 
1512  m  from  the  TM08  simulation.  Although  similar  to 
Fig.  2  in  Terwey  et  al.  (2013),  this  figure  plots  the  evo¬ 
lution  of  the  tangential  and  vertical  velocity  at  the  height 
of  maximum  tangential  winds  in  the  region  where  the 
secondary  eyewall  first  becomes  apparent  in  the  azi- 
muthally  averaged  fields.  As  discussed  in  Terwey  et  al. 
(2013),  Fig.  1  further  evidences  that  the  analyzed  sim¬ 
ulation  exhibits  a  canonical  eyewall  replacement  cycle 
as  inferred  from  observations  (e.g.,  Willoughby  et  al. 
1982;  Houze  et  al.  2007;  Bell  et  al.  2012b;  Sitkowski  et  al. 
2012).  The  RAMS  simulation  is  consistent  also  with  re¬ 
alistic  numerical  integrations  of  hurricanes  that  exhibit 
secondary  eyewalls  (e.g.,  Judt  and  Chen  2010;  Abarca 
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Fig.  1.  Radius-time  plot  of  (a)  azimuthally  averaged  tangential 
velocity  (ms-1)  and  (b)  azimuthally  averaged  vertical  velocity 
(ms-1)  within  200-km  radius  from  vortex  center  from  the  meso- 
scale  model  simulation.  Panels  are  shown  at  1512-m  height,  where 
the  strongest  azimuthally  averaged  tangential  winds  are  found. 
Dashed  lines  highlight  the  three  times  during  the  secondary  eye- 
wall  formation  process  used  often  in  the  subsequent  analysis. 

and  Corbosiero  2011;  Zhou  and  Wang  2011;  Zhou  et  al. 
2011;  Rozoff  et  al.  2012).  Specifically,  at  a  height  of 
1512  m,  the  primary  eyewall  is  centered  around  34-km 
radius  and  its  tangential  wind  intensifies  up  to  hour  15 
(Fig.  la).  Subsequently,  these  tangential  winds  weaken 
monotonically  from  80.9  m  s_1  at  hour  15  to  26.5  m  s_1  at 
hour  42.  As  the  strengthening  and  weakening  occurs  in 
the  primary  eyewall,  the  tangential  winds  outside  the 
primary  eyewall  undergo  a  progressive  amplification  on 
a  relatively  broad  scale  that  extends  beyond  the  domain 
shown  (>200  km).  This  amplification  accounts  for  the 
radial  expansion  of  the  swirling  wind  field  reported  pre¬ 
viously  in  SEF  studies  (e.g.,  Huang  et  al.  2012;  Rozoff 
et  al.  2012).  Superimposed  on  this  broadening  of  the 
tangential  wind,  a  more  localized  tangential  wind  maxi¬ 
mum  emerges  first  around  hour  19  at  85-km  radius  and 
then  subsequently  at  later  times  and  at  other  model 
heights  [discussed  further  in  section  5  and  in  Terwey  et  al. 
(2013)].  This  localized  wind  maximum  continues  to  in¬ 
tensify,  and  exceeds  the  interior  tangential  wind  maxi¬ 
mum  at  about  hour  25  (Fig.  2  of  TM08).  The  new  wind 
maximum  further  intensifies  and  contracts  slowly  inward, 
thereby  completing  an  eyewall  replacement  cycle. 

Figure  lb  shows  the  evolution  of  the  azimuthally  av¬ 
eraged  vertical  velocity  field.  The  strongest  vertical  ve¬ 
locities  attain  values  somewhat  larger  than  1.5 ms-1  and 


are  located  in  the  regions  of  strongest  tangential  winds 
(Fig.  la)  associated  with  the  primary  and  secondary 
eyewalls.  As  the  tangential  wind  velocity  of  the  primary 
eyewall  weakens  (Fig.  la),  so  does  the  corresponding 
vertical  velocity  (Fig.  lb).  In  Fig.  lb,  the  emergence  of  the 
secondary  eyewall  is  reflected  in  the  establishment  of 
a  maximum  of  azimuthal-mean  vertical  velocity  that  oc¬ 
curs  at  nearly  the  same  time  as  the  development  of  the 
localized  tangential  wind  maximum  (Fig.  la).  However, 
near  the  85-km  radius  the  initial  mean  secondary  updraft 
maximum  can  be  recognized  as  early  as  hour  15 — about 
4h  before  the  emergence  of  the  secondary  tangential 
wind  maximum  at  1512-m  height.  In  addition  to  showing 
the  demise  of  the  primary  eyewall  and  the  emergence  of 
the  secondary  eyewall,  Fig.  lb  shows  the  so-called 
moat — a  region  of  pronounced  subsidence  and  greatly 
diminished  radar  reflectivity  (Houze  et  al.  2007)  that  is 
clearly  discernable  between  the  two  eyewalls.  The  sub¬ 
siding  mean  motion  attains  magnitudes  greater  than 
0.5  m  s-1  at  hour  18  near  the  48-km  radius. 

4.  Secondary  eyewall  formation  kinematics 

Figure  2  shows  the  height-radius  structure  of  the  az¬ 
imuthally  averaged  radial  velocity  (u),  tangential  ve¬ 
locity  (v) ,  and  vertical  velocity  (w)  within  200  km  of  the 
circulation  center  during  the  SEF  period.  This  region 
encompasses  both  the  primary  and  secondary  eyewalls. 
The  panels  highlight  characteristic  features  of  the  azi¬ 
muthally  averaged  horizontal  and  vertical  velocity  fields 
that  are  found  during  the  24-h  period  in  which  high- 
time-resolution  output  is  available. 

Figures  2a-c  show  that  strong  inflow  is  confined  to  the 
lowest  km,  with  the  largest  values  at  the  lowest  model 
level  (z  =  147.64  m).  Initially,  the  strongest  values  of 
inflow  are  located  just  outside  the  RMW  of  the  primary 
eyewall.  This  inflow  exceeds  25  m  s-1  at  30-km  radius  at 
hour  4  and  20 ms-1  at  35-km  radius  at  hour  18.  Besides 
the  maximum  inflow  associated  with  the  primary  eyewall 
(>15ms_1),  another  inflow  maximum  progressively 
emerges  in  association  with  the  SEF  shown  in  Fig.  1.  The 
second  inflow  maximum  is  centered  at  about  80-km  ra¬ 
dius  at  18  h  (Fig.  2b)  and  gradually  becomes  the  strongest 
in  the  domain,  exceeding  20 ms-1  at  hour  23  (Fig.  2c). 

Figures  2a-c  show  that  the  main  outflow  resides  at 
upper  levels  (between  11-  and  15-km  height)  and  radially 
exterior  to  the  eyewall  regions.  At  early  times  (Fig.  2a), 
when  the  storm  is  characterized  by  one  eyewall  (Fig.  1), 
the  strongest  outflow  is  more  vertically  confined  than  at 
later  times  (Fig.  2c),  when  the  nascent  secondary  eye- 
wall  projects  into  the  azimuthally  averaged  flow  (Fig.  1). 
In  addition  to  the  upper-level  outflow,  Figs.  2a-c  show 
also  secondary  outflow  regions  located  directly  above 
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Fig.  2.  The  2-h  and  azimuthal  average  of  (a)-(c)  radial,  (d)-(f)  azimuthal,  and  (g)-(i)  vertical  wind  components  (ms-1),  centered  at 
(left)  hour  4,  (middle)  hour  18,  and  (right)  hour  23  from  the  mesoscale  model  simulation  of  TM08.  Contours  are  plotted  every  5  (0.5)  m  s  - 1 
for  the  azimuthal  and  radial  (vertical)  wind  velocities,  except  for  negative  contours  in  the  radial  wind,  which  are  plotted  every  2  ms-1. 
Dashed  lines  indicate  positive  (negative)  values  in  the  radial  (vertical)  wind  component.  The  zero  contour  is  not  plotted.  Vertical  arrows 
denote  the  approximate  radius  of  SEF  at  hours  18  and  23,  as  defined  by  the  secondary  maximum  in  vertical  velocity  in  (h),(i). 


the  inflow  maxima.  As  the  primary  eyewall  weakens,  the 
lower-tropospheric  outflow  associated  with  the  primary 
eyewall  decreases  in  magnitude  from  about  10ms_1  at 
hour  4  to  about  8  m  s-1  at  hour  18  and  to  less  than  5ms-1 
at  hour  23.  In  the  meantime,  the  lower-tropospheric  out¬ 
flow  associated  with  the  SE  increases  in  magnitude  with 
time,  and  by  hour  18  this  outflow  exceeds  4ms-1  (Fig.  2b). 

Figure  2d  shows  that  the  tangential  velocity  structure 
at  hour  4  bears  a  strong  resemblance  to  a  major  tropical 
cyclone  as  observed  with  high-resolution  dropsonde  data 
(e.g.,  Montgomery  et  al.  2006;  Bell  and  Montgomery 
2008).  The  tangential  wind  attains  its  maximum  value 
near  the  top  of  the  boundary  layer  (cf.  Fig.  la)  and 
gradually  weakens  with  height  and  radius.  The  tan¬ 
gential  wind  structure  of  Fig.  2d  evolves  to  a  configu¬ 
ration  where  a  secondary  maximum  is  discernable 
between  70-  and  90-km  radius  at  hour  23  (Fig.  2f).  This 
structure  change  is  associated  with  a  progressive  ra¬ 
dial  broadening  of  the  tangential  wind  field  and  the 
subsequent  development  of  a  secondary  tangential 
wind  maximum  (Fig.  1).  The  broadening  of  the  tangen¬ 
tial  wind  field  occurs  throughout  the  troposphere  and  is 


illustrated  by  inspecting  the  30  and  40  m  s-1  contours  and 
observing  their  progressive  movement  upward  and  out¬ 
wards  with  time  (Figs.  2d-f).  Another  notable  feature  is 
that  the  development  of  the  secondary  tangential  wind 
maximum  occurs  near  the  top  of  the  inflow  layer  (illus¬ 
trated  by  examining  the  evolution  of  the  60  m  s-1  contour 
in  Fig.  2f). 

Figure  2g  shows  (w)  at  hour  4.  At  this  time  the  vertical 
motion  is  confined  mostly  to  the  primary  eyewall  in 
a  relatively  narrow  radial  region  centered  at  about 
30-km  radius  and  sloping  outward  to  about  45-km  radius 
at  9-km  height.  At  later  times  the  vertical  velocity  struc¬ 
ture  of  Fig.  2g  changes  qualitatively  to  a  configuration 
with  a  clear  secondary  maximum  outside  the  primary 
eyewall  (Fig.  2i).  As  the  secondary  eyewall  forms,  a  sec¬ 
ondary  maximum  of  vertical  velocity  emerges  (Fig.  2h) 
and  strengthens  (Fig.  2i)  above  the  region  of  strongest 
tangential  (Figs.  2e,f)  and  radial  (Figs.  2b, c)  winds.  At  early 
times,  the  projection  onto  the  azimuthally  averaged  vertical 
velocity  in  the  nascent  secondary  eyewall  occupies 
a  broader  radial  region  and  possesses  a  weaker  magnitude 
than  the  primary  eyewall.  The  difference  in  the  mean 
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vertical  velocity  structure  between  the  two  eye  walls  may  be 
explained  by  the  fact  that  the  convective  activity  in  the 
secondary  eyewall  occupies  a  broader  radial  region  and  is 
not  contiguous  in  azimuth  (see  Fig.  4  of  TM08). 

During  the  time  when  the  primary  eyewall  tangential 
winds  are  strengthening,  Fig.  2g  shows  that  the  azi- 
muthally  averaged  vertical  velocity  exhibits  two  distinct 
maxima:  one  near  the  top  of  the  inflow  layer  and  another 
one  in  the  middle  troposphere  (centered  at  about  6.5-km 
height).  As  the  primary  eyewall  weakens,  the  maxi¬ 
mum  near  the  top  of  the  inflow  layer  disappears. 
However,  as  the  secondary  eyewall  emerges,  the  two 
maxima  in  the  emerging  eyewall  are  evident  in  Figs.  2h 
and  2i.  The  midtropospheric  updraft  maximum  has 
been  related  to  the  aggregate  radial  gradient  of  the 
buoyancy  generation  rate  resulting  from  the  latent  heat 
release  in  the  deep  cumulus  convection  in  the  central 
region  of  the  vortex  (Hendricks  et  al.  2004;  Montgomery 
et  al.  2006,  2009;  Bui  et  al.  2009).  The  lower  updraft 
maximum  has  been  identified  with  boundary  layer  in¬ 
ertial  effects  in  a  radial  region  where  boundary  layer  in¬ 
flow  terminates  (Smith  et  al.  2009).  The  next  section 
presents  dynamical  evidence  that  the  lower-level  updraft 
maximum  is  indeed  associated  with  inertial  effects,  as  is 
the  tangential  wind  acceleration  near  the  top  of  the  BL 
and  the  emergence  of  the  SE  itself. 

5.  Secondary  eyewall  formation  dynamics 

The  objective  of  this  section  is  to  obtain  a  more 
complete  understanding  of  the  thermodynamical  and 
dynamical  processes  contributing  to  SEF.  We  begin  the 
analysis  with  an  examination  of  the  structure  and  evo¬ 
lution  of  the  azimuthally  averaged  surfaces  of  equivalent 
potential  temperature  6e  and  M.  This  will  furnish  a  zero- 
order  view  of  the  thermodynamical  pathway  to  SEF  in  the 
RAMS  model.  The  illustrated  evolution  of  the  low-level  M 
surfaces,  together  with  the  analysis  in  sections  3  and  4,  then 
motivates  an  analysis  of  departures  of  gradient  wind  bal¬ 
ance  and  a  quantification  of  the  various  flow  contributions 
to  the  tangential  wind  tendency  during  the  SEF  process. 

a.  Evolution  of  0e  and  M  surfaces 

Figure  3  presents  radius-height  plots  of  azimuthal  and 
2-h  averages  of  6e ,  [calculated  using  Eq.  (43)  in  Bolton 
(1980)]  and  M  centered  at  hours  4,  18,  and  23,  as  the 
secondary  eyewall  forms.  These  quantities  are  particu¬ 
larly  informative  since  they  are  approximately  con¬ 
served  following  moist  air  parcels  in  the  free  atmosphere 
assuming  the  azimuthal  pressure  torque  and  small-scale 
mixing  processes  are  negligible. 

Figure  3a  affirms  broadly  the  expected  pattern  of  6e 
within  a  tropical  cyclone  (Hawkins  and  Imbembo  1976; 


Radius  [km] 

Fig.  3.  The  2-h  and  azimuthal  average  of  absolute  angular  mo¬ 
mentum  (rrds1  X  106)  and  equivalent  potential  temperature  (K) 
within  200  km  of  vortex  center  from  the  mesoscale  model  simulation. 
Plots  are  centered  at  (a)  hour  4,  (b)  hour  18,  and  (c)  hour  23.  Contours 
are  plotted  every  1  m2  s_1  X  106  and  5  K.  See  text  for  details. 

Bogner  et  al.  2000),  with  a  markedly  different  structure 
of  the  6e  isopleths  in  the  inner  core  and  radially  outward 
of  it.  In  the  radial  region  where  SEF  takes  place,  and 
particularly  at  radii  greater  than  about  75  km,  the  6e 
surfaces  are  approximately  horizontal.  In  that  radial 
region  there  is  a  decrease  of  6e  with  height  from  the 
surface  up  to  about  3-km  height  (indicating  that,  given 
lifting  and  saturation,  the  lower  troposphere  is  con- 
vectively  unstable  in  the  azimuthal  mean).  Above  3-km 
height,  in  that  same  radial  region,  there  is  an  increase  of 
6e  with  height  (implying  that,  given  lifting  to  saturation, 
the  lower  troposphere  is,  in  the  azimuthal  mean,  con¬ 
vective  unstable).  Closer  to  the  center,  in  the  vicinity  of 
the  primary  eyewall,  the  0e  isopleths  change  their  pat¬ 
tern  by  increasing  their  magnitude  radially  inward  and 
exhibiting  a  structure  weakly  slanted  from  vertical 
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below  about  6  km.  Above  6-km  height,  the  0e  isopleths 
lean  upward  and  outward.  Since  6e  is  essentially  con¬ 
served  above  the  boundary  layer,  the  upward  and  out¬ 
ward  pattern  of  the  6e  isopleths  results  from  the  upward 
and  outward  flow  of  the  primary  eye  wall  (Fig.  2a).  Ra¬ 
dially  inward  of  the  primary  eyewall,  and  as  a  result 
of  strong  subsidence  in  the  eye,  0e  isopleths  exhibit  a 
radius-height  structure  similar  to  that  found  radially 
outside  of  the  eyewall.  As  the  secondary  eyewall  forms, 
the  vertical  motions  are  increasingly  projecting  into  the 
azimuthal  average  (Figs.  2h,i)  and  the  azimuthally  aver¬ 
aged  6e  isopleths  start  to  lean  upward  and  outward 
between  70-  and  140-km  radius. 

The  structure  and  evolution  of  M  isopleths  offer 
dynamical  insight  into  the  SEF  process.  Above  the 
boundary  layer  where  there  is  an  absence  of  strong  dif¬ 
fusion,  the  inward  migration  of  the  M  isopleths  indicates 
an  amplification  of  the  swirling  wind  field  as  discussed 
in  the  introduction.4  A  comparison  of  Figs.  3a-c  evi¬ 
dences  that  the  largest  inward  displacement  of  M  iso¬ 
pleths  occurs  near  the  top  of  the  inflow  layer  (Figs.  2a-c). 
The  top  of  the  strong  inflow  is  where  the  strongest  tan¬ 
gential  wind  velocities  are  observed  and  where  the  sec¬ 
ondary  tangential  wind  maximum  emerges  (Figs.  2d-f). 
Near  the  ocean  surface,  the  M  surfaces  tilt  strongly  in¬ 
ward  with  height.  This  inward  tilt  implies  that  the  tan¬ 
gential  wind  increases  with  height  and  that  M  is  being 
diffused  into  the  ocean  surface  by  nonconservative 
processes.  Near  the  surface,  the  tangential  wind  is  sub¬ 
gradient  relative  to  the  gradient  wind  near  the  top  of  the 
boundary  layer.  Above  the  innermost  point  of  each  M 
surface,  the  M  surfaces  tilt  outward  with  height.  The 
outward  tilt  is  consistent  with  the  ascent  of  air  out  of 
the  boundary  layer  and  its  flaring  out  to  large  radius  in 
the  upper  troposphere  (Fig.  2).  The  fact  that  the  sec¬ 
ondary  eyewall  tangential  maximum  develops  near  the 
top  of  the  BL  (Fig.  2d-f)  suggests  plausibly  that  the 
spinup  of  the  secondary  eyewall  may  be  intimately 
linked  with  the  underlying  boundary  layer  dynamics,  as 
initially  proposed  by  Huang  et  al.  (2012). 

Finally,  in  the  primary  eyewall  (Figs.  2g— i),  the  iso¬ 
pleths  of  6e  and  M  are  nearly  congruent  throughout  the 
troposphere,  with  the  0e  surfaces  sloping  radially  out¬ 
ward.  This  configuration  indicates  that  the  circulation 
tends  to  be  conditionally  symmetrically  neutral 
(Emanuel  1986).  However,  in  the  region  of  SEF,  the  6e 
and  M  isopleths  are  not  congruent  and  intersect  one 
another,  with  the  6e  surfaces  being  approximately 


4  Assuming,  as  discussed  in  the  introduction  and  shown  further 
below,  that  asymmetric  eddy  processes  are  small  and  can  be  ig¬ 
nored  to  a  first  approximation. 


horizontal.  The  evolution  of  the  6e  and  M  surfaces 
demonstrate  that  the  secondary  eyewall  as  simulated 
here  does  not  form  via  a  pathway  of  moist-neutral 
convective  dynamics  (cf.  Nong  and  Emanuel  2003) 
and  points  to  an  intensification  process  (Nguyen  et  al. 
2008;  Montgomery  et  al.  2009;  Smith  et  al.  2009)  that 
is  not  captured  by  the  conceptual  and  modeling 
framework  of  moist-neutral  dynamics  as  described  in 
Emanuel  (1986,  1994,  1995,  1997,  2004). 

b.  Departures  of  balance  in  the  boundary  layer 

To  further  understand  the  boundary  layer  dynamics 
that  contribute  to  the  formation  and  spinup  of  the  sec¬ 
ondary  eyewall,  we  consider  first  an  analysis  of  the 
agradient  force  followed  by  an  analysis  of  the  various 
flow  contributions  to  the  spinup  of  the  tangential  wind  in 
the  boundary  layer  and  above.  The  agradient  force  per 
unit  mass  is  defined  as  the  difference  between  the  local 
radial  pressure  gradient  and  the  sum  of  the  centrifugal 
and  Coriolis  forces  in  the  radial  momentum  equation 
(Smith  et  al.  2009  and  Huang  et  al.  2012): 

Agradient  Force  =  +  fr  +  f(v)  >  (3) 

where  the  symbols  are  as  defined  earlier.  The  agradient 
force  is  a  diagnostic  of  the  departure  of  axisymmetric 
gradient  wind  balance  in  the  radial  coordinate  direction. 
A  positive  agradient  force  implies  an  outward-directed 
acceleration  and  occurs  when  the  tangential  wind  is 
supergradient.  A  negative  agradient  force  implies  an 
inward-directed  acceleration  and  occurs  when  the  flow  is 
subgradient.  Supergradient  flow  indicates  that  the  tan¬ 
gential  wind  is  larger  than  necessary  to  satisfy  gradient 
wind  balance,  while  subgradient  flow  indicates  that  the 
tangential  wind  is  smaller  than  necessary  to  satisfy  gra¬ 
dient  wind  balance. 

Figure  4  shows  radius-height  plots  of  agradient-force 
isopleths  centered  at  4, 18,  and  23  h.  At  all  times  shown  the 
figure  indicates  that  departures  of  gradient  wind  balance 
occur  predominantly  around  the  declining  primary  eyewall 
and  the  emerging  secondary  eyewall  (Fig.  2).  Supergradient 
winds  attain  their  maximum  values  near  the  top  of  the  in¬ 
flow  layer  (Figs.  2a-c)  and  are  collocated  with  the  strongest 
tangential  winds  (Figs.  2d-f)  and  the  largest  inward  dis¬ 
placement  of  the  M  isopleths  (Figs.  3a-c).  As  the  SE  forms, 
the  agradient  force  progressively  increases  in  magnitude.5 
Specifically,  at  hour  4  and  80-km  radius,  the  largest 


5  Consideration  of  a  generalized  Ekman  imbalance  will  render 
consistent  results  with  those  presented  here,  since  the  radial 
frictional  force  in  the  inflow  layer  is  directed  radially  outward. 
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Fig.  4.  The  2-h  and  azimuthal  average  of  agradient  force 
(ms-1h-1)  within  200  km  of  vortex  center  from  the  mesoscale  sim¬ 
ulation.  Plots  are  centered  at  (a)  hour  4,  (b)  hour  18,  and  (c)  hour  23. 
Contours  are  plotted  every  20m s  '  h1.  See  text  for  details. 


agradient  acceleration  is  about  16  m  s_1  h_1,  and  at  hour 
23,  it  surpasses  40ms_1h_1.  Thus,  subgradient  flow  is 
found  radially  outwards  of  the  regions  with  largest 
values  of  supergradient  flow  and  at  the  lowest  level  of 
the  domain  (where  friction  is  the  largest). 

The  results  presented  in  Fig.  4  support  the  findings  and 
interpretations  of  Fluang  et  al.  (2012)  by  confirming  that 
the  supergradient  flow  increases  progressively  within  and 
near  the  top  of  the  boundary  layer  as  the  secondary  wind 
maximum  emerges.  To  further  test  the  proposed  boundary 
layer  control  hypothesis,  we  carry  out  next  a  diagnosis  of 
the  various  flow  contributions  to  the  tangential  velocity 
tendency  equation  during  the  development  of  the  simu¬ 
lated  tangential  wind  maximum. 

c.  Diagnosis  of  tangential  velocity  tendency  equation 

In  cylindrical  coordinates  based  at  the  center  of  cir¬ 
culation,  the  local  time  rate  of  change  of  azimuthally 
averaged  tangential  velocity  can  be  written  as 


d_W) 

dt 


-<“)<£,}  -  WO 


(4) 


where  f a  =  f  +  (1  lr)d(r(v))/dr  is  the  azimuthal-mean 
vertical  component  of  absolute  vorticity.  The  rest  of  the 
symbols  are  as  defined  previously.  On  the  right-hand 
side  of  Eq.  (4),  the  first  two  terms  represent  the  mean 
radial  flux  of  mean  absolute  vorticity  and  the  mean 
vertical  advection  of  (v);  the  next  two  terms  represent 
their  eddy  counterparts — namely,  the  eddy  radial  vor¬ 
ticity  flux  and  the  eddy  vertical  advection  of  eddy  tan¬ 
gential  momentum.  The  fifth  term  represents  the 
azimuthal  perturbation  pressure  gradient  force  per  unit 
mass  and  the  last  term  ((F))  represents  the  combined 
effects  of  surface  friction  and  subgrid-scale  diffusion. 

Figure  5  shows  the  terms  of  Eq.  (4)  for  a  calculation 
centered  on  hour  23,  which  is  the  time  just  prior  to  the 
establishment  of  the  secondary  tangential  wind  maxi¬ 
mum  near  the  top  of  the  boundary  layer  (cf.  Fig.  1).  The 
term  not  shown  is  the  azimuthal  average  of  the  azimuthal 
perturbation  pressure  gradient  force  per  unit  mass,  which 
is  much  smaller  than  the  other  terms  because  of  the 
small  variation  in  the  density  with  the  azimuth.  The 
figure  shows  the  combined  contribution  of  the  two  mean 
terms  (Fig.  5a),  the  combined  eddy  terms  (Fig.  5b),  the 
local  mean  tangential  velocity  tendency  (Fig.  5c),  and 
the  residual  (Fig.  5d).  The  residual  is  calculated  solving 
Eq.  (4)  for  (F)  6 

Figure  5  summarizes  the  contributions  to  the  averaged 
tangential  wind  tendency  for  both  the  spindown  of  the 
primary  eyewall  and  the  spinup  of  the  secondary  eye- 
wall.  The  present  discussion  will  focus  on  the  spinup  of 
the  secondary  eyewall  by  highlighting  the  various  con¬ 
tributions  to  the  tendency  in  the  radial  region  of  the 
SEF,  located  roughly  between  about  75-  and  150-km 
radius.  In  the  lowest  part  of  the  domain  between  about 
75  and  150  km  (Fig.  2f),  mean  tangential  wind  velocity 
change  achieves  values  of  about  lms-1h-1  (Fig.  5c). 
This  acceleration  is  the  result  of  the  substantially  larger 
and  opposing  contributions  of  different  terms  in  the 
tangential  velocity  budget.  The  terms  with  the  largest 
magnitudes  are  the  combined  mean  terms  (Fig.  5a)  and 


6  In  the  RAMS  output,  the  friction  term  is  not  readily  available 
and  so  the  residual  serves  as  an  estimate  of  the  frictional  force  per 
unit  mass.  This  term  contains  not  only  the  effect  of  frictional  forces, 
but  also  inaccuracies  associated  with  the  interpolation  from 
the  staggered  Arakawa  C  grid  (in  the  Cartesian  coordinates  of  the 
RAMS  simulation)  to  the  cylindrical  coordinate  grid  used  in  the 
diagnostic  analyses. 
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Fig.  5.  The  2-h  average  of  the  terms  in  the  tangential  momentum 
equation  centered  around  hour  23  from  the  mesoscale  simulation. 
The  terms  shown  are  (a)  the  combined  contribution  of  the  mean 
radial  vorticity  flux  and  the  mean  vertical  advection  of  mean  tan¬ 
gential  momentum,  (b)  the  combined  contribution  of  the  pertur¬ 
bation  radial  eddy  vorticity  flux  and  the  perturbation  vertical 
advection  of  perturbation  tangential  momentum,  (c)  the  local  time 
change  of  tangential  velocity,  and  (d)  the  residual,  estimated  as  the 
difference  between  (c)  and  the  sum  of  (a)  and  (b).  For  clarity, 
contours  are  shown  every  2ms”1h”1  in  (a),(b),  0.5ms_1h_1  in 
(c),  and  5ms-1h-1  in  (d).  The  zero  contour  is  not  shown. 

the  residual  (Fig.  5d).  The  residual  is  associated  with 
nonconservative  effects  owing  to  surface  friction  in  the 
azimuthal  direction,  and  is  maximized  in  the  frictionally 
boundary  layer  near  the  surface. 

The  combined  mean  terms  (Fig.  5a)  contribute  to 
a  positive  tendency  mainly  below  2  km  and  between 
about  60-  and  150-km  radius.  By  definition,  the  mean 
advection  terms  are  the  sum  of  both  the  radial  flux  of 
mean  absolute  vorticity  and  the  mean  vertical  advection 
of  (v).  In  the  boundary  layer,  the  largest  values  of  these 
two  terms  arise  from  the  radial  influx  of  mean  absolute 
vorticity  and  this  influx  exceeds  30ms_1h_1  (not 
shown).  The  mean  vertical  advection  of  (v)  attains 
values  of  more  than  —  10ms_1h_1  (not  shown)  and 
serves  to  loft  the  enhanced  tangential  wind  generated  in 
the  boundary  layer  to  the  vortex  aloft.  The  mean  vertical 
advection  term  thus  serves  as  a  negative  tendency 


“sink”  in  the  boundary  layer  and  a  positive  tendency 
“source”  in  the  vortex  aloft  between  1-  and  2-km  height. 

In  the  SEF  region,  the  combined  effect  of  the  eddy 
terms  to  the  tangential  wind  tendency  is  small  overall 
and  negative  in  the  boundary  layer  (Fig.  5b).  Examining 
the  individual  components  of  the  eddy  terms  (not  shown), 
we  find  that  the  eddy  radial  flux  of  eddy  absolute  vorticity 
in  the  region  of  SEF  is  largely  negative.  On  the  other 
hand,  the  vertical  eddy  advection  of  eddy  tangential 
momentum  accounts  for  a  negative  tangential  wind 
tendency  below  2  km  and  a  positive  tendency  above  that 
height.  The  latter  thus  contributes  to  a  spinup  of  the 
tangential  wind  above  the  boundary  layer. 

The  foregoing  results  demonstrate  that  the  actual  SEF 
tangential  velocity  tendency  in  the  SEF  region  occurs 
within  but  near  the  top  of  the  boundary  layer  even  in  the 
presence  of  a  strong  spindown  tendency  associated  with 
surface  friction.  The  results  show  that  the  dominant 
contribution  to  the  spinup  of  the  mean  tangential  wind  is 
associated  with  the  mean  radial  influx  of  absolute  vor¬ 
ticity  in  the  boundary  layer  of  the  vortex.  The  results 
indicate  also  that  in  the  region  of  SEF  the  main  effect  of 
the  eddy  vertical  advection  is  to  redistribute  (v) ,  serving 
to  spin  down  (v)  in  the  lowest  1  or  2  km  and  to  spin  up  (v) 
in  a  layer  above  that.  Together,  these  findings  confirm  the 
operation  of  the  boundary  layer  spinup  mechanism  in  the 
region  of  SEF  as  foreshadowed  in  the  introduction. 

Although  the  foregoing  results  provide  strong  support 
for  the  new  boundary  layer  control  hypothesis,  a  linger¬ 
ing  question  remains  as  to  whether  there  is  a  clear  causal 
link  between  the  boundary  layer  dynamics  and  SEF. 

d.  Causality  and  the  arrow  of  time 

Figure  6  shows  the  time  evolution  of  the  tangential 
wind  velocity,  the  horizontal  divergence  and  the  agra- 
dient  force  within  the  boundary  layer  at  the  model-level 
height  of  457  m.7  For  each  quantity  presented,  both 
eyewalls  are  recognizable  as  independent  maxima.  The 
primary  eye  wall  is  centered  between  25-  and  30-km  ra¬ 
dius  and  the  SE,  which  becomes  increasingly  clear  with 
time,  exhibits  a  radial  contraction  as  described  previously 
in  Fig.  1. 


7  Taking  a  boundary  layer  average  of  the  lowest  962-m  results  in 
a  figure  consistent  with  Fig.  6,  particularly  in  its  space-time  struc¬ 
ture  (not  shown).  The  conclusions  presented  in  this  manuscript 
thus  hold  true  with  both  approaches  (i.e.,  single-level  or  layer- 
averaged  views).  The  main  difference  between  the  two  approaches 
is  that  the  vertical  average  comprises  larger  subgradient  wind 
magnitudes,  which  is  to  be  expected  because  of  the  inclusion  of  the 
lowest  levels  that  possess  the  largest  deficit  from  the  gradient  flow 
at  the  boundary  layer  top. 
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Fig.  6.  Radius-time  plot  of  azimuthal  averages  of  (a)  tangential 
wind  velocity  (ms-1),  (b)  horizontal  divergence  (s  _1  X  103),  and 
(c)  agradient  force  (m  s-1  h-1)  at  457-m  height  from  the  mesoscale 
simulation.  Color  contours  are  plotted  every  20  ms  _1h“\  s_1  X  103, 
and  5 ms-1,  respectively.  Red  horizontal  arrows  and  hour  numbers 
indicate  the  earliest  time  that  a  secondary  maximum  of  the  plotted 
variable  is  discernible  in  the  azimuthal  average.  [For  reference, 
hour  23,  the  time  of  the  earliest  tangential  velocity  secondary 
maximum,  at  1512  m  (see  Fig.  1)  is  indicated.] 

Figure  6b  shows  the  time  evolution  of  horizontal 
divergence  within  the  boundary  layer.  Horizontal 
convergence  exhibits  a  secondary  maximum  in  the 
boundary  layer  beginning  around  14  h  at  65-km  radius 
and  beginning  around  15  h  at  85-km  radius.  Recalling 
the  chronological  evolution  presented  in  Figs,  la,  2d-f, 
and  3a,  the  emergence  of  the  convergence  maximum 
precedes  the  emergence  of  the  secondary  tangential 
wind  maximum  at  any  height.  The  secondary  conver¬ 
gence  maximum  at  85-km  radius  appears  simultaneously 
with  the  secondary  maximum  in  vertical  velocity 
documented  in  Fig.  1.  We  will  now  show  that  although 
the  secondary  maximum  in  convergence  and  vertical 
velocity  occur  nearly  simultaneously,  their  emergence 
is  preceded  by  the  maximum  in  the  agradient  force 
field. 

Figure  6c  presents  the  evolution  of  the  agradient  wind 
force  field  within  the  boundary  layer.  It  shows  that  the 


agradient  force  field  develops  a  secondary  maximum  at 
hour  13  at  65-km  radius  and  at  hour  14  at  85-km  radius. 
This  secondary  maximum  precedes  by  1  h  the  emergence 
of  a  secondary  maximum  in  radial  convergence  and  ver¬ 
tical  velocity.  Moreover,  this  secondary  maximum  in  the 
agradient  force  field  precedes  by  10  h  the  emergence  of 
the  secondary  tangential  wind  maximum  (near  the  top  of 
the  boundary  layer).  This  chain  of  events  suggests 
strongly  that  the  positive  agradient  force,  which  quan¬ 
tifies  the  effective  outward  radial  force  in  the  boundary 
layer,  is  the  cause  for  the  increased  radial  convergence 
and  vertical  motion  on  the  system  scale  in  the  region  of 
SEF.  This  sequence  of  events  supports  the  idea  that  di¬ 
agnostics  of  departures  of  gradient  wind  balance  in  the 
boundary  layer  may  prove  useful  as  a  SEF  forecasting 
tool  (Huang  et  al.  2012).  We  look  now  for  another  test  of 
the  boundary  layer  control  hypothesis  using  an  in¬ 
dependent  modeling  platform. 

6.  Further  tests  of  the  boundary  layer  control 
hypothesis  using  a  simple  slab  model 

In  the  previous  sections  it  was  shown  that  the  emer¬ 
gence  of  the  secondary  wind  maximum  occurs  near  the 
top  of  the  boundary  layer  and  follows  a  chain  of  events 
initiated  by  a  substantial  radial  expansion  of  the  tan¬ 
gential  wind  field  (associated  with  the  first  spinup 
mechanism)  and  the  development  of  supergradient 
winds  in  the  boundary  layer  (associated  with  the  second 
spinup  mechanism).  The  analysis  of  the  SEF  kinematics 
and  dynamics  in  the  RAMS  simulation  was  thus  shown 
to  be  consistent  with  the  new  SEF  model  developed  by 
Huang  et  al.  (2012),  who  proposed  that  the  underlying 
boundary  layer  dynamics  of  a  maturing  hurricane  vortex 
are  an  important  controlling  element  in  SEF.  This  idea  is 
quite  different  from  other  SEF  hypotheses  summarized 
in  the  introduction.  As  with  any  new  scientific  idea,  it 
should  be  tested  from  different  perspectives.  The  au¬ 
thors  of  this  manuscript  are  further  testing  the  boundary 
layer  control  hypothesis  using  other  full-physics  hurri¬ 
cane  simulations  as  well  as  observational  datasets.  The 
results  of  these  analyses  will  be  reported  in  due  course. 
Here,  we  offer  a  complimentary  approach  for  testing  the 
proposed  relationship  between  boundary  layer  dynam¬ 
ics  and  SEF. 

The  test  is  a  relatively  simple  one:  If  SEF  is  largely  a 
boundary  layer  response  to  a  radially  growing  tangential 
wind  field,  one  would  expect  that  a  boundary  layer 
model  would  respond  suitably  in  generating  a  secondary 
wind  maximum  using  a  fixed,  simulated  tangential  wind 
field  above  the  boundary  layer  in  the  full-physics  model, 
together  with  reasonable  estimates  of  the  drag  coefficient, 
the  boundary  layer  depth,  and  horizontal  eddy  diffusivity. 
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Vsiab-BL  and  Uslab_BL 

- VHAMS  and  UHAMS  at  z  =  1138  m 

- VHAMS  and  UHAMS  at  z  =  787  m 

- VHAMS  and  Urams  at  z  =  457  m 

- -Vrams  and  Urams  at  z  =  148  m 


Fig.  7.  Tangential  and  radial  wind  field  profiles  of  the  mesoscale  simulation  after  azimuthal  and  2-h  averages  at  the  z  =  1138-,  787-,  457-, 
and  148-m  model  levels  (dashed  curves)  and  slab  boundary  layer  simulations  (instantaneous  fields,  solid  black  curves),  (a)  Wind  profiles  at 
hour  13  of  the  mesoscale  simulation  including  the  tangential  wind  profile  at  z  =  1512  m  that  is  used  as  the  initial  condition  for  the  slab 
model  and  forcing  (solid  curve;  see  text  for  details),  (b)  Mesoscale  model  profiles  centered  at  hour  18  and  the  results  of  the  slab  model 
(solid  curve)  after  5  h  of  integration,  (c)  Mesoscale  model  profiles  centered  at  hour  23  and  the  results  of  the  slab  model  (solid  curve)  at  10  h 
of  integration,  when  the  slab  model  is  reaching  equilibrium. 


To  execute  such  a  test  we  integrate  the  slab  boundary 
layer  model  described  in  section  2b,  using  as  initial 
condition  the  RAMS  azimuthally  averaged  tangential 
wind  profile  at  the  top  of  the  boundary  layer.  The  wind 
field  is  taken  at  hour  13,  prior  to  the  time  that  SE  is 
evident  in  any  of  the  azimuthally  averaged  fields  of 
Fig.  1.  We  then  compare  the  results  of  the  RAMS  and 
slab  boundary  layer  models.  To  illustrate  how  the  slab 
model  typically  achieves  the  computed  equilibrium,  radial 
profiles  of  the  tangential  and  radial  velocity  are  shown  at 
the  slab  model  times  of  5  and  10  h.  The  solid  curves  in 
Figs.  7b  and  7c  show  the  temporal  evolution  of  the  slab 
boundary  layer  model,  which  reaches  an  equilibrium  at 
about  hour  11  of  the  integration.  For  reference,  these 
outputs  are  compared  with  the  5-h  output-interval  data 
from  the  RAMS  simulation  centered  at  hours  13, 18,  and 
23  (dashed  lines;  Figs.  7a-c,  respectively)  using  an  azi¬ 
muthal  and  2-h  average  of  the  tangential  and  radial  wind 
velocities  at  z  —  1138-,  787-,  457-,  and  148-m  model  levels, 
respectively.  Also  shown  in  Fig.  7a  is  the  initially  zero 
radial  inflow8  and  correspondingly  averaged  tangential 
wind  from  the  RAMS  simulation  at  z  =  1.512  km,  both  of 
which  were  used  to  integrate  the  slab  boundary  layer 
model  forward  in  time  (dark  solid  curves). 

Figure  7  shows  that  the  slab  boundary  layer  develops 
a  secondary  maximum  in  the  tangential  wind  field.  It 
shows  also  that  such  a  maximum  is  not  qualitatively 
unlike  that  found  in  the  RAMS  simulation.  Figure  7b 
demonstrates  that  the  secondary  tangential  wind  maximum 


8  Integrations  of  the  slab  boundary  layer  model  initialized  with 
the  radial  inflow  from  the  RAMS  simulation  exhibit  a  very  similar 
evolution  to  those  initialized  with  zero  radial  inflow  after  the  first 
hour  or  so  when  a  rapid  inertial  adjustment  of  the  inflow  occurs. 


in  both  models  emerges  clearly  after  only  5h  of  nu¬ 
merical  integration  beyond  13  h  (the  initial  condition  of 
the  slab  boundary  layer  integration).  At  the  time  shown 
in  Fig.  7b,  the  secondary  tangential  wind  maximum  in 
the  slab  model  is  centered  at  105  km  and  reaches  a  value 
of  50 ms-1.  The  RAMS  secondary  tangential  wind 
maximum  is  centered  at  about  85  km  with  a  magnitude 
of  55  m  s-1  at  z  =  1.138  km.  As  in  the  typical  evolution  of 
eye  wall  replacement  cycles  (e.g.,  Willoughby  et  al.  1982; 
Willoughby  and  Black  1996),  the  secondary  tangential 
wind  maximum  in  both  models  becomes  more  pro¬ 
nounced  with  time  as  this  feature  contracts  radially.  At 
the  time  of  Fig.  7c,  the  slab  boundary  layer  tangential 
wind  maximum  attains  a  maximum  value  of  58 ms-1, 
centered  at  80-km  radius.  The  RAMS  tangential  wind 
maximum  reaches  a  value  of  60 ms-1  at  80-km  radius. 
Figures  7b  and  7c  show  that,  as  in  the  case  of  the  tan¬ 
gential  wind,  the  radial  velocity  of  the  slab  boundary 
layer  model  is  qualitatively  consistent  with  the  corre¬ 
sponding  radial  profiles  of  the  radial  velocity  found  in 
the  RAMS  simulation.  Both  models  exhibit  radial  inflow 
with  two  maxima,  and  each  maximum  resides  just  out¬ 
side  the  corresponding  tangential  wind  maximum. 

The  development  of  a  secondary  wind  maximum  is 
found  to  be  a  robust  feature  in  the  slab  boundary  layer 
model.  This  has  been  verified  with  a  variety  of  sensitivity 
experiments  using  initial  conditions  and  tangential-wind 
forcing  taken  from  different  vertical  levels  and  times 
around  hour  13  of  the  RAMS  simulation  (not  shown).  In 
all  experiments  conducted,  the  secondary  tangential  wind 
maximum  emerges  consistently  with  comparable  magni¬ 
tudes  and  radial  location.  These  sensitivity  experiments 
demonstrate  that  given  a  prescribed  tangential  wind  field 
at  the  top  of  the  boundary  layer  and  reasonable  estimates 
of  the  surface  drag  coefficient,  radial  eddy  diffusivity, 
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and  boundary  layer  depth,  axisymmetric  boundary  layer 
dynamics  alone  are  capable  of  developing  a  secondary 
tangential  wind  maximum  and  a  concomitant  region  of 
focused  convergence. 

Although  the  boundary  layer  tangential  and  radial 
velocities  in  the  two  models  exhibit  qualitative  consis¬ 
tency,  there  are  quantitative  differences  among  them.  In 
the  slab  model,  for  example,  the  tangential  wind  maxi¬ 
mum  and  its  corresponding  radial  wind  structure  are 
somewhat  more  pronounced  than  in  the  RAMS  model. 
Also,  the  tangential  wind  maximum  in  the  slab  model 
has  a  more  localized  secondary  tangential  wind  maxi¬ 
mum  than  that  of  the  RAMS  model.  However,  given  the 
high  degree  of  idealization  invoked  in  the  slab  boundary 
layer  model  and  its  intrinsic  limitations  (e.g.,  Smith  and 
Vogl  2008),  quantitative  differences  between  the  models 
are  to  be  expected. 

The  imposition  of  a  fixed  tangential  wind  profile  prior 
to  the  emergence  of  the  simulated  tangential  wind 
maximum  is  not  necessarily  the  most  accurate  way  to 
emulate  the  changing  wind  field  above  and  within  the 
boundary  layer  during  the  SEF  process.  However,  this 
methodology  has  shown  that  a  secondary  tangential 
wind  maximum  can  emerge  solely  as  a  result  of  bound¬ 
ary  layer  dynamics.  Furthermore  these  experiments 
produced  an  unexpected  result;  namely,  the  simulated 
secondary  wind  maximum  contracts  inwards,  as  sec¬ 
ondary  eyewalls  in  mesoscale  simulations  and  in  nature 
do.  Specifically,  the  inward  contraction  of  the  secondary 
wind  maximum  in  the  slab  model  occurs  without  special 
constraints  on  the  underlying  surface  heat  fluxes  or 
special  inertial  stability  considerations  on  the  balanced 
secondary  circulation  above  the  boundary  layer  (cf. 
Emanuel  1994, 1997;  Nong  and  Emanuel  2003;  Schubert 
and  Hack  1982) — properties  that  are  irrelevant  to  the 
slab  boundary  layer  model  used  here.  These  findings  are 
interesting  because  they  point  to  a  hitherto  unrecognized 
role  of  unbalanced  boundary  layer  dynamics  in  the  in¬ 
ward  contraction  of  the  outer  eyewall  under  an  imposed 
wind  field  at  the  boundary  layer  top. 

7.  Summary 

In  this  paper  we  have  sought  to  obtain  an  improved 
understanding  of  secondary  eyewall  formation  by  a  care¬ 
ful  analysis  of  the  phenomenon  in  two  different  modeling 
frameworks.  The  first  of  them  is  a  high-resolution,  full- 
physics,  three-dimensional  mesoscale  model  in  which  we 
study  an  integration  that  undergoes  a  canonical  eyewall 
replacement  cycle.  For  this  simulation,  it  is  shown  that 
secondary  eyewall  formation  occurs  in  a  region  of  the 
vortex  where  azimuthal  averages  of  equivalent  potential 
temperature  and  absolute  angular  momentum  surfaces 


are  far  from  being  congruent.  This  result  questions  the 
applicability  of  moist-neutral  viewpoints  and  related 
mathematical  formulations  thereto  for  studying  this 
process  of  tropical  cyclone  intensity  change. 

The  analysis  of  the  mesoscale  simulation  fully  supports 
the  recent  findings  of  Huang  et  al.  (2012).  An  inves¬ 
tigation  of  the  kinematics  of  secondary  eyewall  formation 
shows  that  the  establishment  of  the  secondary  eyewall  is 
the  culmination  of  a  sequence  of  structural  changes  in 
the  tropical  cyclone.  The  sequence  begins  with  a  radial 
broadening  of  the  tangential  wind  field,  followed  by  the 
emergence  of  persistent  convective  activity  in  the  sec¬ 
ondary  eyewall  formation’s  radial  region,  and  subse¬ 
quently  followed  by  the  emergence  of  the  secondary 
wind  maximum  near  the  top  of  the  boundary  layer.  An 
analysis  of  the  dynamics  of  the  secondary  eyewall  for¬ 
mation  process  confirms  that  the  broadening  of  the 
tangential  wind  field  above  the  boundary  layer  and  the 
associated  increase  in  the  radial  pressure  gradient 
induces  a  boundary  layer  response  that  includes  a  rela¬ 
tively  large  radial  inflow  radially  outside  the  primary 
eyewall.  This  strong  inflow  results  in  an  outward-directed 
agradient  force  that  causes  sustained  horizontal  conver¬ 
gence.  This  convergence,  in  the  conditionally  unstable 
region  where  the  secondary  eyewall  forms,  leads  to  sus¬ 
tained  deep  convective  activity.  The  strong  inflow  there 
results  also  in  a  strong  mean  influx  of  absolute  vorticity 
that  surpasses  the  depletion  of  absolute  angular  mo¬ 
mentum  by  surface  friction.  When  the  radial  inflow  has 
attained  sufficient  strength,  the  net  effect  of  the  vorticity 
influx  and  frictional  loss  in  the  boundary  layer  is  the 
generation  of  a  secondary  tangential  wind  maximum 
within  the  boundary  layer  outside  the  primary  eyewall. 

The  second  modeling  framework  is  an  axisymmetric, 
nonlinear,  time-dependent,  slab  boundary  layer  model 
with  radial  diffusion.  When  this  model  is  forced  with  the 
azimuthally  averaged  pressure  field  (associated  with  the 
gradient  wind)  at  the  top  of  the  boundary  layer  from 
the  mesoscale  simulation  at  times  well  before  the  vortex 
exhibits  a  secondary  wind  maximum,  the  model  generates 
a  secondary  tangential  wind  maximum.  The  secondary 
wind  maximum  is  qualitatively  comparable  in  both 
magnitude  and  time  evolution  to  the  maximum  predicted 
by  the  mesoscale  simulation.  These  results  demonstrate 
that  axisymmetric  boundary  layer  dynamics  alone  are 
capable  of  developing  secondary  wind  maxima.  The  re¬ 
sults  suggest  also  that  unbalanced  boundary  layer  dy¬ 
namics  may  make  an  important  quantitative  contribution 
to  the  contraction  of  the  tangential  wind  maximum. 
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